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O : Abstract 

^ • We observe Bose-Einstein correlations in 7r° pairs produced in Z° hadronic decays using 

Qh! the data sample collected by the OPAL detector at LEP 1 from 1991 to 1995. Using a 

r-| ! static Gaussian picture for the pion emitter source, we obtain the chaoticity parameter A = 
0.55 ± 0.10 ± 0.10 and the source radius R = (0.59 ± 0.08 ± 0.05) fm. According to 

■ the JETSET and HERWIG Monte Carlo models, the Bose-Einstein correlations in our data 

■ sample largely connect 7r°s originating from the decays of different hadrons. Prompt pions 
' formed at string break-ups or cluster decays only form a small fraction of the sample. 



(Submitted to Physics Letters B) 



The OPAL Collaboration 



G. Abbiendi^, C. Ainsley^, P.F. Akesson^, G. Alexander^^, J. AUison^^, P. Amaral^, 
G. Anagnostou\ KJ. Anderson^, S. Arcelli^, S. Asai^^, D. Axen^'^, G. Azuelos^^-", 
1. Bailey2^ E. Barberio^f , R.J. Barlowl^ R.J. Batley^ R Bechtle^^ T. Behnke^-', 
K.W. Belpo, RJ. Bell\ G. Bella^^, A. Bellerive^ G. Benelli^ S. Bethke^^, Q. BiebeP^, 
I.J. Bloodworth\ O. Boeriu^°, P. Bock^\ D. Bonacorsi^, M. Boutemeur^\ S. Braibant^, 

L. Brigliadori^, R.M. Brown^°, K. Buesser^^, H.J. Burckhart^, S. Campana^, 
R.K. Carnegie^ B. Caron^s, A.A. Carter^^^ J.R. Carter^ C.Y. Chang^^ D.G. Charltoni'\ 
A. CsillingS'^ M. Cuffiani^, S. Dado^i, A. De Roeck^ E.A. De Wolf^ K. Desch^^, 

B. Dienes^o, M. Donkers^ J. Dubbert3\ E. Duchovni^^, G. Duckeck3\ I.R Duerdothl^ 
E. Elfgrenl^ E. Etzion^^, R Fabbri^, L. Feld^°, R Fe^■ari^ R Riedler3\ I. Heck^o, M. Rord^ 

A. Rrey^ A. Rurtjes^ R Gagnon^^ J.W. Gary^, G. Gaycken^^ C. Geich-Gimbel^, 
G. Giacomelli^, R Giacomelli^, M. Giunta^ J. Goldberg^i, E. Gross^^ J. Grunhaus^^ 
M. Gruwe^ RO. Gunther^ A. Gupta^ C. Hajdu^^, M. Hamann^^, G.G. Hanson\ 
K. Harder^^, A. HareP\ M. Harin-Dirac'^, M. Hauschild^, CM. Hawkes\ R. Hawkings^, 

R.J. Hemingway^ C. Hensel^^, G. Herten^o, R.D. Heuer^^, J.C. HilP, K. Hoffman^ 
R.J. Horner^, D. Horvattf^''^, P. Igo-Kemenes^^, K. Ishii^'^, H. Jeremie^^, P. Jovanovic\ 
T.R. Junk^, N. Kanaya^^, J. Kanzaki^^, G. Karapetian^^ D. Karlen*', K. Kawagoe^^ 
T. Kawamoto^^ R.K. Keeler^^, R.G. Kellogg^^ B.W. Kennedy^o, D.H. Kim^^ K. Klein^^'*, 
A. Klier2\ S. Kluth'^^^ T. Kobayashi^^, M. Kobel'\ S. Komamiya^^, L. Kormos^e, 
T. Kramer2^ T. Kress^, R Krieger^ ', J. von Kroghi\ D. Krop^^ K. Kruger^ T. KuhP^, 
M. Kupper^^, G.D. Lafferty^^, H. Landsman^i, D. Lanske^'', J.G. Layter^, A. Leins^\ 
D. Lellouch^^, J. Letts", L. Levinson^^, J. Lillichi", S.L. Lloyd^^ RK. Loebinger^s, J. Lu^^ 
J. Ludwig^°, A. Macpherson^^'*, W. Mader^, S. Marcellini^, A.J. Martin^^, G. Masetti^, 

T. Mashimo^^, R Mattig™, W.J. McDonald^^, J. McKenna^^ T.J. McMahon\ 
R.A. McPherson^^, R. Meijers*, W. Menges2^ RS. Merritt^, H. Mes^'", A. Michelini^, 
S. Mihara^^ G. Mikenberg^^, D.J. Miller^^ S. Moed"^^, W. Mohr^", T. Mori^^, A. Mutter^", 
K. Nagai^^ I. Nakamura^^, H.A. NeaF, R. Nisius^^, S.W. 0'Neale\ A. Oh^ A. Okparai\ 
M.J. Oreglia^, S. Orito^^, C. Pahl■■^^ G. Pasztor^ ^ J.R. Paterl^ G.N. Patrick^o, J.E. Pilcher^ 

J. Pinfold^^ D.E. Plane^ B. Poli^, J. Polok^ O. Pooth^^ M. Przybycieri^ ", A. Quadt^, 
K. RabbertzS'^ C. Rembser^ R RenkeP^, H. Rick^ J.M. Roney^^, S. Rosati^, Y. Rozen^^ 

K. Runge^°, K. Sachs^, T. Saeki^^, E.K.G. Sarkisyan^'^ A.D. Schaile^\ O. Schaile^\ 
P. Scharff-Hansen^, J. Schieck'^^, T. Schorner-Sadenius^, M. Schroder^, M. Schumacher^, 

C. Schwick^ W.G. Scott^o, R. Seuster^^'-^, T.G. Shears^'^, B.C. Shen\ R Sherwood^^ 
G. Siroli^ A. Skuja^^ A.M. Smith^ R. Sobie^^, S. Soldner-Rembold^^''^, R Spano^ 

A. Stahl^, K. Stephens^^, D. Strom^^, R. Strohmer3\ S. Tarem^^ M. Tasevsky^ 
R.J. Taylor^^ R. Teuscher^, M.A. Thomson^, E. Torrence^'^, D. Toya^^, P. Tran^, 
A. Tricoli^, I. Trigger^, Z. Tr6csanyi^°'*^, E. Tsur^^, M.R Turner- Watson^, I. Ueda^'^, 
B. Ujvari3°'% C.R VoUmer^i, R Vannerem^", R. Vertesi=^°, M. Verzocchi^^ H. Yoss^'i, 

J. Vossebeld^'^, D. Waller^ C.R Ward^ D.R. Ward^ RM. Watkins\ A.T. Watson\ 
N.K. Watson\ RS. Wells^ T. Wengler^ N. Wennes^ D. Wetterling^^ G.W. Wilson^^.fc^ 
J.A. Wilson\ G. Wolf^^, T.R. Wyattl^ S. Yamashita^^, D. Zer-Zion^ L. Zivkovic^^ 



^School of Physics and Astronomy, University of Birmingham, Birmingham B 15 2TT, UK 
^Dipartimento di Risica dell' Universita di Bologna and ENRN, 1-40126 Bologna, Italy 



1 



^Physikalisches Institut, Universitat Bonn, D-53115 Bonn, Germany 

"^Department of Physics, University of California, Riverside CA 92521, USA 
^Cavendish Laboratory, Cambridge CBS OHE, UK 

^Ottawa-Carleton Institute for Physics, Department of Physics, Carleton University, Ot- 
tawa, Ontario KIS 5B6, Canada 

^CERN, European Organisation for Nuclear Research, CH-1211 Geneva 23, Switzerland 
^Enrico Fermi Institute and Department of Physics, University of Chicago, Chicago IL 

60637, USA 

^°Fakultat fiir Physik, Albert-Ludwigs-Universitat Freiburg, D-79104 Freiburg, Germany 
^^Physikalisches Institut, Universitat Heidelberg, D-69120 Heidelberg, Germany 
^^Indiana University, Department of Physics, Bloomington IN 47405, USA 

Queen Mary and Westfield College, University of London, London El 4NS, UK 
^^Technische Hochschule Aachen, III Physikalisches Institut, Sommerfeldstrasse 26-28, 
D-52056 Aachen, Germany 

i^University College London, London WCIE 6BT, UK 

^^Department of Physics, Schuster Laboratory, The University, Manchester M13 9PL, UK 
^^Department of Physics, University of Maryland, College Park, MD 20742, USA 
^^Laboratoire de Physique Nucleaire, Universite de Montreal, Montreal, Quebec H3C 3J7, 
Canada 

^^University of Oregon, Department of Physics, Eugene OR 97403, USA 
20CLRC Rutherford Appleton Laboratory, Chilton, Didcot, Oxfordshire OXll OQX, UK 
^^Department of Physics, Technion-Israel Institute of Technology, Haifa 32000, Israel 
^^Department of Physics and Astronomy, Tel Aviv University, Tel Aviv 69978, Israel 
^^International Centre for Elementary Particle Physics and Department of Physics, Univer- 
sity of Tokyo, Tokyo 1 13-0033, and Kobe University, Kobe 657-8501, Japan 

Particle Physics Department, Weizmann Institute of Science, Rehovot 76100, Israel 
^^Universitat Hamburg/DESY, Institut fiir Experimentalphysik, Notkestrasse 85, D-22607 
Hamburg, Germany 

26University of Victoria, Department of Physics, P O Box 3055, Victoria BC V8W 3P6, 
Canada 

^^University of British Columbia, Department of Physics, Vancouver BC V6T IZl, Canada 
^^University of Alberta, Department of Physics, Edmonton AB T6G 2J1, Canada 
^^Research Institute for Particle and Nuclear Physics, H-1525 Budapest, P O Box 49, Hun- 
gary 

^"Institute of Nuclear Research, H-4001 Debrecen, P O Box 51, Hungary 
^^Ludwig-Maximilians-UniversitatMiinchen, Sektion Physik, Am Coulombwall 1, D-85748 

Garching, Germany 

^^Max-Planck-Institute fiir Physik, Fohringer Ring 6, D-80805 Miinchen, Germany 
^^Yale University, Department of Physics, New Haven, CT 06520, USA 



" and at TRIUMF, Vancouver, Canada V6T 2A3 

and Royal Society University Research Fellow 

and Institute of Nuclear Research, Debrecen, Hungary 

and Heisenberg Fellow 
^ and Department of Experimental Physics, Lajos Kossuth University, Debrecen, Hungary 
^ and MPI Miinchen 



2 



^ and Research Institute for Particle and Nuclear Physics, Budapest, Hungary 
now at University of Liverpool, Dept of Physics, Liverpool L69 3BX, U.K. 
^ and CERN, BP Div, 1211 Geneva 23 

^ now at University of Nijmegen, HEFIN, NL-6525 ED Nijmegen,The Netherlands, on 
NWO/NATO Fellowship B 64-29 

^ now at University of Kansas, Dept of Physics and Astronomy, Lawrence, KS 66045, 
U.S.A. 

' now at University of Toronto, Dept of Physics, Toronto, Canada 

current address Bergische Universitat, Wuppertal, Germany 
" and University of Mining and Metallurgy, Cracow, Poland 
° now at University of California, San Diego, U.S.A. 

P now at Physics Dept Southern Methodist University, Dallas, TX 75275, U.S.A. 
' now at IPHE Universite de Lausanne, CH-1015 Lausanne, Switzerland 
^ now at lEKP Universitat Karlsruhe, Germany 

* now at Universitaire Instelling Antwerpen, Physics Department, B-2610 Antwerpen, Bel- 
gium 

* now at RWTH Aachen, Germany 



1 Introduction 



The Bose-Einstein correlations (BEC) effect has a quantum-mechanical origin. It arises 
from the requirement to symmetrise the wave function of a system of two or more identical 
bosons. It was introduced into particle reactions leading to multi-hadron final states as 
the GGLP effect |1| in the study of the 7r+7r+ and tt^tt^ systems. The distributions of 
the opening angle between the momenta in pairs of like- sign pions were shifted towards 
smaller values compared to the corresponding distributions for unlike-sign pairs. A related 
effect was exploited earlier in astronomy [2 J to measure the radii of stars. 

In high energy physics, for example e+e^ collisions at LEP, a quantitative understand- 
ing of the BEC effect allows tests of the parton fragmentation and hadronisation models. 
This would in turn help in achieving a more precise measurement of the W boson mass 
and better knowledge of several Standard Model (SM) observables [3 1. The fragmentation 
models presently used are those of strings and clusters implemented, respectively, in the 
JETSET [i4i and HERWIG [5| Monte Carlo generators. 

Numerous studies of BEC in pairs of identical bosons already exist, see for example 0. 
Due to the experimental difficulties in photon and 7r° reconstruction, only very few stud- 
ies lITl exist for BEC in 7r° pairs, even though they offer the advantage of being free of final 
state Coulomb corrections. 

The string model predicts a larger BEC strength or chaoticity and a smaller effective 
radius of the emitting source for 7r° pairs compared to vr^ pairs while the cluster fragmenta- 
tion model predicts the same source strength and size [ 8 , 9 1 . However, neither model of pri- 
mary hadron production has a mechanism to allow BEC between 7r°s produced in different 
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strong decays. The string model prediction is a consequence of electric charge conserva- 
tion in the local area where the string breaks up. Similar expectations can be derived if the 
probabilities in the string break-up mechanism are interpreted as the squares of quantum 
mechanical amplitudes [8, 10 1. A small difference between tt^ pairs and 7r° pairs is also 
expected from a pure quantum statistical approach to Bose-Einstein symmetry ifTTl . In ad- 
dition, based on isospin invariance, suggestions exist on how to relate BEC in the pion-pair 
systems i.e. tt^tt^, tt^tt^, and tt+tt^ and how to extend it to t^^-k^. IfT^l . The L3 collabo- 
ration has recently reported |7| that the radius of the neutral-pion source may be smaller 
than that of charged pions, Rt,±t,± — Rt^OttO = (0.150 ± 0.075(stat.) ± 0.068(syst.)) fm, in 
qualitative agreement with the string fragmentation prediction. 

This paper presents a study of BEC in 7r° pairs using the full hadronic event sample 
collected at centre-of-mass energies at and near the Z° peak by the OPAL detector at LEP 
from 1991 to 1995. This corresponds to about four million hadronic Z° decays. A highly 
pure sample of 7r° mesons is reconstructed using the lead-glass electromagnetic calorimeter. 
The correlation function is obtained after accounting for purity and resonant background. 
It is parametrised with a static picture of a Gaussian emitting source HIEI. 



2 Selection of hadronic decays 

A full description of the OPAL detector can be found in llT3l . The sub-detectors relevant to 
the present analysis are the central tracking detector and the electromagnetic calorimeter. 
The central tracking detector consists of a silicon micro-vertex detector, close to the beam 
pipe, and three drift chamber devices: the vertex detector, a large jet chamber and surround- 
ing 2r-chambers^ In combination, the three drift chambers sitting inside a solenoidal mag- 
netic field of 0.435 T yield a momentum resolution of o-pjpt ~ a/0.022 _j_ (0.0015 ■ pt)^ 
for I cos(^^)| < 0.7, where pt (in GeV) is the transverse momentum with respect to the 
beam axis. The electromagnetic calorimeter detects and measures the energies and posi- 
tions of electrons, positrons and photons for energies above 0.1 GeV. It is a total absorbing 
calorimeter, and is mounted between the coil and the iron yoke of the magnet. It consists 
of 1 1704 lead-glass blocks arranged in three large assemblies (the barrel that surrounds the 
magnet coil, and two endcaps) which together cover 98% of the solid angle. The intrinsic 
energy resolution is cr^/i? ~ 5%/ ^/E, where E is the electromagnetic energy in GeV. 

Standard OPAL selection criteria are applied to tracks and electromagnetic clusters (X^. 
Tracks are required to have at least 20 measured points in the jet chamber, a measured mo- 
mentum greater than 0.1 GeV, an impact parameter \do\ in the r — (p plane smaller than 
2 cm, a z position at the point of closest approach to the origin in the r — (p plane within 25 
cm of the interaction point, and a measured polar angle with respect to the beam axis greater 
than 20°. Electromagnetic clusters are required to have an energy greater than 0.1 GeV if 
they are in the barrel part of the detector (i.e. | cos^^ | < 0.82) or greater than 0.3 GeV if 
they are in the endcap parts. Hadronic Z° decays are selected by requiring for each event 

^The OPAL coordinate system is defined so that the z axis is in the direction of the electron beam, the x 
axis points towards the centre of the LEP ring, and 9 and </> are the polar and azimuthal angles, defined relative 
to the +z- and +a;-axes, respectively. In cylindrical polar coordinates, the radial coordinate is denoted r. 
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more than 7 measured tracks, a visible energy larger than 60 GeV and an angle larger than 
25° and smaller than 155° between the calculated event thrust 111311 axis and the beam axis. 
The visible energy is the energy sum of all detected tracks, electromagnetic clusters not 
associated to tracks and electromagnetic clusters associated to tracks after correcting for 
double counting. A sample of 3. 1 million Z° hadronic decays is selected for which the total 
background, consisting mainly of r pairs, is less than 1% and is neglected throughout the 
analysis. 

Detector effects and detection efficiencies for the spectra of 7r° pairs are evaluated using 
eight million Monte Carlo hadronic Z° decays. Events are generated using the JETSET 7.4 
program, tuned to reproduce the global features of hadronic events as measured with the 
OPAL detector (14], with the BEC effect explicitly switched off. Samples generated with 
the HERWIG 5.9 program without the BEC effect are used for comparison. The generated 
events were passed through a full simulation of the OPAL detector [ 16 1 and were analysed 
using the same reconstruction and selection programs as were applied to the data. 



3 Reconstruction of tt" mesons 

For the selected event sample, neutral pions are reconstructed from photon pairs. Pho- 
ton reconstruction is performed in the barrel part of the electromagnetic calorimeter where 
both the photon reconstruction efficiency and the energy resolution are good. The pro- 
cedure of [17 1 which resolves photon candidates in measured electromagnetic clusters is 
used. It employs a parametrisation of the expected lateral energy distribution of electro- 
magnetic showers. It is optimised to resolve as many photon candidates as possible from 
the overlapping energy deposits in the electromagnetic calorimeter in a dense environment 
of hadronic jets. The purity of the photon candidate sample is further increased using a 
likelihood-type function [ 17 1 that associates to each photon candidate a weight w for being 
a true photon. Photon candidates with higher w are more likely to be true photons. 

All possible pairs of photon candidates are then considered. Each pair was assigned 
a probability P for both candidates being correctly reconstructed as photons. This prob- 
ability is simply the product of the w -weights associated with the two candidates. The 
combinatorial background consists of a mixture of three components: (i) wrong pairing of 
two correctly reconstructed photons, (ii) pairing of two fake photons and (iii) pairing of 
one correctly reconstructed photon with a fake one. Choosing only photon pairs with high 
values of P leaves combinatorial background mostly from component (i). 

The 7r° reconstruction efficiency and purity are illustrated in Figure [T]for different cuts 
on P. The efficiency is defined as the ratio of the number of correctly reconstructed tt^s over 
the number of generated 7r°s , and the 7r° purity is defined as the ratio of signal over total 
entries in a photon-pair mass window between 100 and 170 MeV. 
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Figure 1 : The tt" reconstruction efficiency (top) and the purity (bottom) for different cuts 
on the weight P = wi x wj of the ij photon pair. The purity and efficiency are estimated 
from the JETSET Monte Carlo. The corresponding statistical errors are smaller than 1%. 
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4 Selection of tt" Pairs 



The average number of 7r°s produced in Z° decays has been measured ifTSll to be 9.76±0.26, 
which is reproduced by our Monte Carlo simulations. This leads to about 45 possible 
7r° pairings per event. Considering only 7r° candidates with P > 0.1 (i.e. 17% efficiency 
and 36% purity), we reconstruct at the detector level 4.7 7r° candidates on average per 
event. This leads to about 8 pairings among which only 1 pair on average is really formed 
by true 7r°s . Here, the detector level means that detector response, geometrical acceptance 
and photon reconstruction efficiency are taken into account. Therefore, the tt" pair sample 
is background dominated and the study of 7r° pair correlations or invariant mass spectra is 
subject to very large background subtraction. Monte Carlo must be used to predict both the 
shape and amount of background to be subtracted, leading to large systematic errors in the 
measurements of the BEC source parameters. 

To avoid this, the 7r° selection criteria are tightened. We select vr^s which have a mo- 
mentum above 1 GeV. This cut reduces the fraction of fake tt^s . In addition, it removes 
7r°s produced by hadronic interactions in the detector material for which the Monte Carlo 
simulation is not adequate. The probability P associated to each 7r° candidate is required 
to be greater than 0.6. In the case where a photon can be combined in more than one 
pair, only the pair with the highest probability is considered as a 7r° candidate. Among the 
events with four or more reconstructed photon candidates, only those leading to a possible 
7r° pair with four distinct photon candidates are retained for further analysis. Events with 
six or more photon candidates leading to more than two 7r° candidates are rejected. They 
represent about 10% of the retained sample and would increase the sensitivity to unwanted 
resonance signals if they were not rejected. Figure |2l shows the photon pair mass, M2^, 
for the selected events. The average purity of the 7r° sample is 79% in the mass window 
between 100 and 170 MeV. The background is estimated directly from data by a second 
order polynomial fit to the side bands of the peak and by Monte Carlo simulation. The two 
background estimations yield compatible results and the Monte Carlo reproduces correctly 
the data. The superimposed curves are not the result of a fit to the data, but smoothed his- 
tograms of the Monte Carlo expectations for signal and background normalised to the total 
number of selected hadronic Z° decays. 

A clear 7r° pair signal is obtained as shown in Figure |3] where the two values of 
are shown for the retained events. A 7r° pair is considered as a signal candidate if both 
values of are within the mass window between 100 and 170 MeV. The average 7r° pair 
signal purity is 60% and the Monte Carlo simulation describes the data well. Kinematic 
fits were made, constraining the mass of pairs of photon candidates to the vr^ mass, with 
the assumption that the photons come from the primary interaction vertex. Monte Carlo 
studies showed that this gives a 26% improvement in the resolution of the 7r° momentum. 
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Figure 2: Distribution of two-photon invariant mass, M2'y, for selected events which have 
exactly two reconstructed 7r° candidates per event. The smooth curves represent the total 
Monte Carlo expectation (solid line) and the background (dashed line) expectation. The 
curves are normalised to the same number of total selected hadronic Z° decays as in the 
data. The shaded region (100-170 MeV) represents the selected window for the 7r° signal. 
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5 The BEC Function 



The correlation function is defined as the ratio, 

where Q is a Lorentz-invariant variable expressed in terms of the two 7r° four-momenta pi 
and p2 via = — (pi — ^2)^, piQ) = {l/N)dN/dQ is the measured Q distribution of 
the two 7r°s and po{Q) is a reference distribution which should, in principle, contain all the 
correlations included in p{Q) except the BEC. For the measurement of po{Q), we consider 
the two commonly used methods |16| : 



• Event Mixing: Mixed pairs are formed from vr^s belonging to different Z° decay 
events in the data. To remove the ambiguity on how to mix events, we select two-jet 
events having a thrust value T > 0.9, i.e. well defined back-to-back two-jet events. 
The thrust axes of the two events are required to be in the same direction within 
(Acos6' X A0)= (0.05 x 10°). Mixing is then performed by swapping a 7r° from 
one event with a 7r° from another event. To avoid detection efficiency problems 
arising from different detector regions, swapping of two pions is performed only if 
they point to the same region of the electromagnetic barrel detector within (A cos 9 x 
A0)= (0.05 X 10°). With this procedure, we start with two hadronic Z° events each 
having two 7r° candidates and can end up with between zero and four pairs of mixed 
7r° candidates. The Q variable is then calculated for each of the mixed pairs. If 
the contributions from background are removed or suppressed, this method offers 
the advantage of being independent of Monte Carlo simulations, since C (Q) can be 
obtained from data alone. 

• Monte Carlo Reference Sample: The po distribution is constructed from Monte Carlo 
simulation without BEC. The Monte Carlo is assumed to reproduce correctly all 
the other correlations present in the data, mainly those corresponding to energy- 
momentum conservation and those due to known hadron decays. In order to be 
consistent with the first method, the cut T > 0.9 is also applied for both data and 
Monte Carlo. 



In the following, the distributions p{Q) and po(Q) are measured from the same sample of 
selected events. The mixing technique is used as the main analysis method and the Monte 
Carlo reference technique is applied only for comparison. 



6 The Measured BEC Function and Background Contri- 
bution 



The correlation function, C{Q), corresponds experimentally to the average number of 
7r° pairs, corrected for background, in the data sample divided by the corresponding cor- 
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rected average number in the reference sample. Thus, we can write 

^t^. = p(Ql = p'"(Q)-/(Q) ... 
PoiQ) pW)-p^(Q)' ^ 

where and are the measured values, and and Pq are the corresponding corrections 
for background contributions. For both the numerator and denominator, the background 
consists mainly of 7r° pairs in which one or both 7r° candidates are fake. 

The background distributions p^ and pg are obtained from the Monte Carlo information. 
These background distributions can also be obtained from data using a side band fit to the 
projected spectra of the two-dimensional distributions (see FigureEl) in each 400 MeV 
interval of the measured Q variable. The resulting background distributions are correctly 
reproduced by Monte Carlo. However, for the smaller Q intervals as used in this analysis, 
i.e. 100 MeV, the side band fit is subject to large statistical fluctuations, so the Monte Carlo 
distributions have to be used. 

In the region of interest where the BEC effect is observed, Q < 700 MeV, pion 
pairs from particle (resonance) decays could mimic the effect. The relevant decays are: 
K° -> 7r°7r°, fo(980) 7r°7r°, and rj 7r°7r°7r° with branching ratios of 39%, 33% and 
32% respectively. Pion pairs from 7] decay contribute only to the region Q < 315 MeV. 
According to Monte Carlo studies, the number of reconstructed Kg in the 2 7r° channel is 
very small. Furthermore, the hypothesis that each 7r° originates from the primary vertex, as 
used in the kinematic fits (Section 4), does not apply. This is an advantage for this analysis 
since the peak is flattened, making its effect on the Q distribution negligible. The Monte 
Carlo estimates of this particle decay backgrounds are included in the distribution p^{Q), 
adjusting the rate of individual hadrons to the LEP average ifTHI where necessary. 

For our analysis we select 7r° candidates with momentum greater than 1 GeV This is 
dictated by the observation of correlations at small Q even for Monte Carlo events gen- 
erated without any BEC effect. Indeed, as shown in Figure HI a clear BEC-type effect 
is visible in the correlation function obtained from Monte Carlo events without BEC for 
different low cuts on 7r° momentum. Using Monte Carlo information, we find that these 
correlations are mainly caused by tt^s originating from secondary interactions with the de- 
tector material. They would constitute an irreducible background to the BEC effect if low 
momentum vr^s are considered in the analysis. This effect vanishes for 7r° momenta greater 
than 1 GeV. 

We rely on Monte Carlo simulation only to define the appropriate momentum cut (i.e. 
1 GeV) which completely suppresses the effect of soft pions produced in the detector ma- 
terial, rather than relying on its prediction for the exact shape and size of this effect. The 
reason is that, in contrast to charged pions where the measured track information can be 
used to suppress products of secondary interactions in the detector material, the neutral 
pions have to be assumed to originate from the main interaction vertex. Furthermore, with 
this assumption the kinematic fits (Section 4) bias the energy of soft pions emitted in the 
detector material towards larger values since the real opening angle between the photons is 
larger (vertex closer to the calorimeter) than the assumed one. 

With the above selection criteria, the composition of the selected 7r° pair sample is 
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studied using Monte Carlo simulations. According to the string fragmentation model im- 
plemented in JETSET, the selected sample consists of about 97.9% of mixed pion-pairs 
from different hadron decays, 2% of pairs belonging to the decay products of the same 
hadron and only 0.1% prompt pairs from the string break-ups. Similarly, using the clus- 
ter fragmentation model implemented in HERWIG, the selected sample consists of 97% 
of pairs from different hadron decays, 2.3% belonging to the decay products of the same 
hadron and only 0.7% originating directly from cluster decays. It is worth mentioning that 
even if the direct pion pairs from string break-up (JETSET) or cluster decays (HERWIG) 
were all detected and accepted by the analysis procedure, they would be diluted in combi- 
nation with other pions and would constitute only a marginal fraction (< 1% ) of the total 
number of reconstructed 7r° pairs. Thus, our analysis has no sensitivity to direct pion pairs 
originating from string break-up or cluster decay. 



7 Results 

The correlation distribution C (Q) (Eq. @) is parametrised using the Fourier transform of 
the expression for a static sphere of emitters with a Gaussian density (see e.g. fll9il ): 

CiQ) = N{1 + Xexp{-R''Q^)){1 + 5Q + eQ^). (3) 

Here A is the chaoticity of the correlation [which equals zero for a fully coherent (non- 
chaotic) source and one for a chaotic source], R is the radius of the source, and a nor- 
malisation factor. The empirical term, (1 + 5Q + eQ^), accounts for the behaviour of the 
correlation function at high Q due to any remaining long-range correlations. The C{Q) 
distribution for data is shown in Figure|5]as the points with corresponding statistical errors, 
and the smooth curve is the fitted correlation function in the Q range between and 2.5 
GeV. A clear BEC enhancement is observed in the low Q region of the distribution. The 
parameters are determined to be: 

A= 0.55 ±0.10, 

R= (0.59 ± 0.08) fm, 

N= 1.10 ±0.08, 

S= (-0.14 ± 0.05) GeV-\ 

e= (0.07 ± 0.03) GeV-2, 

where the quoted errors are statistical only and the x^/ndf of the fit is 14.7/19. 

The distribution C{Q) obtained for Monte Carlo events generated with no BEC is 
shown as a histogram in the same figure. It shows that there is no residual correlation 
at low Q and indicates that the observed enhancement is present in the data only. The 
dashed-line histogram of Figure |5] represents the correlation function obtained from data 
but before the subtraction, using the Monte Carlo estimates, of pairs from the decay prod- 
ucts of the same hadron, indicating that these contributions have only a minor influence 
on the measured parameters. In addition, the correlation function constructed with back- 
ground 7r° pairs does not show any enhancement at low Q (not shown). Here, background 
7r° pairs are defined as pairs for which one or both of the vr^s are outside the mass window 
100-170 MeV, i.e. these are likely to be fake 7r° candidates. 
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Figure 4: The correlation distribution C{Q) determined for JETSET Monte Carlo events 
(generated without BEC effect) for different cuts on the 7r° momenta, p^. 



13 



The second method, which uses the MC reference sample, yields the following results: 

A =0.50 ±0.10, 
R = (0.46 ± 0.08)fm. 

These results are quoted for comparison only. We choose to quote the results obtained with 
the event mixing method since they are much less dependent on details of the Monte Carlo 
modelling. 

The string model predicts a smaller source radius and a larger chaoticity in the BEC 
effect for 7r° pairs than for n"^ pairs, while the cluster model predicts no difference. These 
predictions hold only for prompt boson pairs produced directly from the string or cluster 
decays. According to our Monte Carlo simulations, we have no sensitivity to these pairs. 



8 Systematic Uncertainties 

Potential sources of systematic error are investigated. In each case the effect on the param- 
eters R and A and their deviations with respect to the standard analysis are estimated. The 
results are summarised in Table[ll 

• Bin width resolution: After the kinematic fits (Section 4), the resolution on the 
invariant mass of two pions, or on the variable Q, is approximately 60 MeV. We have 
chosen a bin width of 100 MeV for the fit to the measured C{Q) distribution. This 
bin width is varied from 100 MeV to 80 MeV and to 120 MeV. 

• Fit range: The low end of the fit range is set to start at Q = 350 MeV (fourth bin) 
The high end of the fit range is changed to stop at Q = 2 GeV 

• Effect of hadron decays: To estimate the effect of the 7r° pairs from the same 
resonance decay on the measured BEC parameters, the estimated contribution is var- 
ied by ±10% which represents the typical error on the measured individual hadron 
rates ifTSl . In order to investigate the dependence of the measured parameters R and 
A on the tt" momentum cut, the analysis is repeated for 7r° momenta larger than 1 .2 
GeV. 

• Analysis procedure: The analysis is repeated for several variations of the selection 
criteria. 

- 1) The 7r° selection mass window is changed from 100-170 MeV to 110-165 
MeV ( increases the tt^ purity by 5%). 

- 2) The probability for 7r° selection is changed from 0.6 to 0.5 (reduces the 
7r° purity by 5%). 

- 3) The thrust value for two-jet events is changed from 0.9 to 0.85 and to 0.92 
(changes the overall event sample size by ±5%). 

- 4) The factor 1 + 6Q + eQ^ is replaced by 1 ± ^Q. 
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Figure 5 : The correlation distribution C(Q) as measured for OPAL data. The smooth curve 
is the fitted correlation function and the dotted histogram is the correlation distribution 
obtained for JETSET Monte Carlo events generated without BEC. The dashed histogram 
represents the measured correlation function before the subtraction of the contributions 
from known hadron decays. 
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- 5) 7r° from different events are mixed if they point to the same region of the 
detector within (A cos 6* x A0) = (0.10 x 15°) instead of (0.05 x 10°). 



The final systematic errors are obtained from quadratically adding the deviations from the 
central value. Thus, 

A = 0.55 ±0.10 ±0.10, 

R = (0.59 ± 0.08 ± 0.05) fm. 



9 Conclusions 



We have observed Bose-Einstein correlations of n pairs produced in hadronic Z decays. 
Assuming a Gaussian shape for the source, we obtain A = 0.55 ± 0.10 ± 0.10 for the 
chaoticity parameter and R = (0.59 ± 0.08 ± 0.05) fm for the radius. In order to con- 
struct a reference sample with the event mixing method, this analysis is restricted to well 
defined back-to-back two-jet events. Furthermore, in order to remove 7r°s not originating 
from the primary interaction vertex the considered momentum phase space is restricted to 
p^o > 1 GeV. The measured value of the source radius is smaller than our former value [I20II . 
R = (1.002±0.016l 

ami) f'^' obtained for charged pions for which the measured track pa- 
rameters allowed access to lower momenta and where the reference sample was constructed 
with unlike-sign pion pairs. However, the value is compatible with the LEP inclusive av- 
erage lEDl . R = (0.74 ± 0.01 ± 0.14) fm, for charged pions. Pions from strong decays 
constitute the dominant part of our sample of reconstructed vr^ pairs. We have no sensi- 
tivity to test the string or cluster model predictions concerning differences between neutral 
and charged pion pairs. We deduce that Bose-Einstein correlations exist between 7r° pairs 
in which each 7r° is a strong decay product of a different hadron. 
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